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ABSTRACT 

We report observations of a sample of high redshift sources (1.8 < z < 4.7), mainly radio-quiet quasars, 
at 350 fj,m using the SHARC bolometer camera at the Caltech Submillimeter Observatory. Nine sources 
were detected (> Aa) and upper limits were obtained for 11 with 350 /im flux density limits (3(t) in the 
range 30 - 125 mJy. Combining published results at other far-infrared and millimeter wavelengths with the 
present data, we are able to estimate the temperature of the dust, finding relatively low values, averaging 
50 K. From the spectral energy distribution, we derive dust masses of a few 10^ H^qqMq and luminosities of 
4-33 X 10^^ H^qqLq (uncorrected for any magnification) implying substantial star formation activity. Thus 
both the temperature and dust masses are not very different from those of local ultraluminous infrared 
galaxies. For this redshift range, the 350 /im observations trace the 60-100 /im rest frame emission and 
are thus directly comparable with IRAS studies of low redshift galaxies. 

Subject headings: infrared: galaxies — quasars: individual (BR12G2— 0725, BRI1335— 0417, 4C41.17, 
H1413-H117, F10214-F472) — infrared: ISM: continuum 
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1. INTRODUCTION 

The study of molecular gas and dust is a significant ob- 
servational tool for probing the physical conditions and 
star formation activity in local galaxies. Recent advances 
in observational techniques at submillimeter and millime- 
ter wavelengths now permit such studies to be made at 
cosmological distances. 

As a result of the IRAS survey many local objects are 
recognized as containing a large mass of dust and gas, such 
that the objects may be more luminous in the far-infrared 
(FIR) than in the optical. The question as to how many 
such objects there may be at cosmological distances, and 
whether they can ac count for the recently discovered FIR 



cosmic background (Puget et al. 1996; Fixsen et al. 1998 
Hauser et al. 1998), is attracting much interest and spawiT 



mg new ins t rument construction and new surveys ( Hughei 
et al. 19981; |Ivison ct al. 19981 ; iKawara et al. 1998| ; pujel 

T^; Lilly et a 



et al. 19991 



arger et al. 199^ ; [Lilly et al. 1999 ). Existing 
submillimeter cameras on ground-based telescopes are not 
yet sensitive and large enough to detect distant objects at 
350 ^m in arbitrary blank fields, e.g., the initial Caltech 
Submillimeter Observatory (CSO) SHARC survey which 
achieved 100 m Jy (1 a) over about 10 square arcminutes 
. However, such cameras can measure the 
objects of known position. A step forward 



dPhiUips 1997| ) 
350 /im flux of 
in the field was the recognition of IRAS F10214+4724 as 
a high redshift object (z — 2.286) by Rowan-Robinson et 
al. (1991). However, it has proved difficult to find many 
such objects to study. On the other hand, quasars have 
sometimes proved to exist in the environs of dusty galaxies 
( Haas et al. 199^ ; Lewis et al. 1998 ; Downcs fc Solomoij 
1998). Omont et al. (1996b) have shown by means of a 



1300 /im survey of radio-quiet quasars that the dust emis- 
sion at high redshifts can be detected in a substantial frac- 
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tion of their project sources. 

In this letter we present measurements at 350 /im to- 
wards a sample of 20 sources with redshifts 1.8 < z < 4.7 
which were selected from different surveys and studies (ref- 
erences are given m Tables I & 2). The sources had 
previously been detected at longer wavelengths, and are 
predominantly from the work of Omont et al. (1996b) 
and Hughes, Dunlop & Rawlings (1997). A wavelength 
of 350 /fm roughly corresponds to the peak flux density of 
highly redshiftcd (z ~ 3) dust emission of objects with 
temperatures of 40 to 60 K. Together with measurements 
at longer wavelengths, it strongly constrains the dust tem- 
perature and, hence, the dust mass and, especially, the 
luminosity of the object. Some of these results were first 
presented by Benford et al. (1998). In this paper, we 
assume Hq — /iioo x lOOkms"^ Mpc~^ and JIq = 1. 



2. OBSERVATIONS AND RESULTS 

The measurements were made during a series of observ- 
ing runs in 1997 February and October and 1998 January 
and April with the 10.4 m Leighton telescope of the CSO 
on the summit of Mauna Kea, Hawaii, during excellent 
weather conditions, with 225 GHz atmospheric opacities of 
< 0.05 (corresponding to an opacity of 1.5 at 350 /xm). 
We used the CSO bolometer camera, SHARC, described 
by Wang et al. (1996) and Hunter, Benford, & Serabyn 
(1996). It consists of a linear 24 element close-packed 
monolithic silicon bolometer array operating at 300 mK. 
During the observations, only 20 channels were opera- 
tional. The pixel size is 5" in the direction of the array and 
10" in the cross direction. The weak continuum sources 
were observed using the pointed observing mode with the 
telescope secondary chopping in azimuth by 30" at a rate 
of 4 Hz. The telescope was also nodded between the on 
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TABLE 1 

Sources detected at 350 and derived properties. 



Source 


z 




R.A. 


Dec. 




Flux Density 


Tdust 


Dust Mass* 


Luminosity* 


Ref. 










(B1950.0) 






(mjy, ±1 (t) 


(K) 


(lOS/iio^o Mo) 


(10l2h-2 Lq) 




BR1202-0725 


4.69 


12 


02 


49.3 


-07 


25 


50 


106±7 


50±7 


4Q+0.9 


14.91°;? 


1,2,3,4 


BRI1335-0417 


4.41 


13 


35 


27.6 


-04 


17 


21 


52±8 


43±6 


•^•"-0.8 




1,2,5 


HMOOOO-263 


4.10 


DO 


00 


49.5 


-26 


20 


01 


134±29 


60t 




20.7+1-g 


2,6 


4C41.17 


3.80 


06 


47 


20.8 


+41 


34 


04 


37±9 


52±6 


1 o+°i 




7,8 


PC2047+0123 


3.80 


20 


47 


50.7 


+01 


23 


56 


80±20 


50t 


r, O + 0.6 

^•■^-0.6 




9,10 


Q1230+1627 


2.74 


12 


30 


39.4 


+16 


27 


26 


104±21 


49±12 


9 c: + 0.4 


S 9+2-0 


2,11 


QOlOO+1300 


2.68 


01 


00 


33.4 


+13 


00 


11 


131±28 


68±5 


1 9+0.5 
1-^-0.3 


23.8tO;2 


12,13 


H1413+117* 


2.54 


14 


13 


20.1 


+11 


43 


38 


293±14 


45±7 


8 9+1-^ 




14,15 


F10214+4724* 


2.28 


10 


21 


31.1 


+47 


24 


23 


383±51 


55±3 




32.8l6;5 


16 



t Because of a lack of sufficient data at other wavelengths, the temperatures of these sources are poorly constrained. The available upper limits 
constrain HMOOOO— 263 to 50 < T(J^gt < 75K, so we adopt 60K; PC2047+0.123 is essentially unconstrained, so we adopt 50K as that is roughly 
the median value for the dust temperature. 

■I-The dust masses and luminosities listed above have not hoen co rrected for any lensing amplific ation. The probable intrinsic value s may be 
found by dividing by A = 7.6 for H 1413+117 (AUoin et al. 1997) and A = 13 for F10214+4724 (Downes, Solomon, & Radford 1995). 



References.— 1. Storrie-Lombardi et al. (1996); [2.] Omont et al. (1996b); [3.] McMahon et al. (1994); [4.] Isaak et al. (1994); [5.] Guil- 
loteau et al. (1997); 6. Schneider, Schmidt & Gunn (1989); 7. Chambers, Miley & van Bruegel (1990); [8.] Hughes, Dunlop & Rawhngs(1997); 
9. Schneider, Schmidt & Gunn (1994); [10.] Ivison (1995); 11. Foltz et al. (1987); 12. Steidel & Sargent (1991); [13.] Guilloteau et al. (1999); 
14. Hazard et al. (1984); [15.] Barvainis, Coleman & Antonucci (1992); [16.] Rowan-Robinson et al. (1993) 
Note: the references in brackets refer to the observations shown in Fig. 2. 



and the off beams at a rate of ~0.1Hz. The point-source 
sensitivity of SHARC at 350 /xm is ^ lJy/-\/Hz and the 
beam size is ^ 9" FWHJVI. All measurements were made 
at 350/im with the exception of H1413-fll7 which was 
also observed at 450/im. 

Pointing was checked regularly on nearby strong galac- 
tic sources which also served as secondary calibrators, and 
was found to be stable with a typical accuracy of <J 3". 
The planets JVIars, Saturn and Uranus served as primary 
flux calibrators. The absolute calibration was found to 
be accurate to within 20%. Repeated observations of 
H 1413 -t- 117 and F 10214 4724 confirmed a relative flux 
accuracy of ^ 20%. The data were reduced using the 
CSO BADRS software package. Typical sensitivities (Icr) 
of ~ 20 mJy were achieved, after ~ 2500 s of on-source 
integration time. 

Nine sources were detected at levels of 4 cr and above, 
as outlined in Table |l|. Included are the z > 4 quasars 
BR 1202-0725, BRI 1335-0417 an d HM 0000-263. Ex- 
cept f or the Clove rleaf (H 1413-1-117; Barvainis, AntonucciJ 
& Co eman 1992), the present measurements are the first 
reported detections for high redshift quasars at 350 /im. 
IVIany of the sources were measured three or more times 
providing both consistency checks and improvements in 
the accuracy of the flux densities. The two strongest 
sources, H 1413-^117 and IRAS F10214-I-4724, were often 
measured before starting the long (~2-3 hours) integra- 
tions on the weaker sources. 

As an illustration of the data quality, Figure 1 
shows the 350 im\ CSO-SHARC measurement towards 
BR 1202— 0725 a.t z — 4.69. This measurement corre- 
sponds to a total of 4 hours integration on source. The 



source is centered at offset zero. The other channels pro- 
vide a measure of the neighboring blank sky emission and 
a reference for the quality of the detection. 

11 sources with redshifts between 1.8 and 4.5 were not 
detected at 350 ^m, with flux density upper limits at the 
3(7 level of 30-125mJy. Table 2 lists their names, red- 
shifts, 350 fj,m flux density measurements with ±1 a errors, 
and a 3 (T upper limit to their luminosities (see below). 




-50 50 

Offset (arcseconds) 

Fig. 1 - The 350 /im flux density measured in the bolometers 
of the SHARC array towards BR1202-0725 at z=4.69. Off- 
sets are given in arcsec with respect to the reference channel 
number 10. The source is centered at offset zero. Emission is 
also seen in the two neighboring pixels because the bolometers 
sample the diffraction pattern of the telescope with a Nyquist 
sampling. 



3 



TABLE 2 
Sources with upper limits at 350 fj.M 



Source 


z 


S350 


i^FlR [ocr) 


Ref. 






(mjy, ±1(t) 


iW^''h,^fo L©) 




BR2237-0607 


4.56 


5 ± 15 


< 6 


1 


BRI0952-0115 


4.43 


8± 22 


< 8 


1,2 


PSS0248+1802 


4.43 


-75 ±22 


< 9 


3 


BR1117-1329 


3.96 


27± 13 


< 4 


1 


Q0302-0019 


3.28 


34 ±21 


< 5 


4 


Q0636+680 


3.18 


-123 ±38 


< 9 


5 


Q2231-0015 


3.01 


-65 ± 24 


< 6 


4 


MG0414+0534 


2.64 


-24 ±35 


< 8 


6 


Q0050-2523 


2.16 


69 ±42 


< 9 


4 


Q0842+3431 


2.13 


16± 10 


< 2 


1 


Q0838+3555 


1.78 


39± 19 


< 4 


1 



References. — 1. Omont et al. (1996b); 2. Guilloteau et al. 
(1999); 3. Kennefick et al. (1995); 4. Hewett, Foltz & Chaffee 
(1995); 5. Sargent, Steidel & Boksenberg (1989); 6. Barvainis et 
al. (1998) 



3. DISCUSSION 

Figure |^ displays the spectral energy distributions of 
the six sources detected at 350 for which fluxes at two 
or more other wavelengths are available from the litera- 
ture. In the following, we will first comment on individual 
sources and then discuss the physical properties of the ob- 
jects. 

The two radio-quiet, z > 4 quasars, BR 1202-0725 
and BRI 1335—0417, are exceptional objects with large 
masses of gas (~ 10^^ M0) which have been detected in 
CO by Omont et al. (1996a) and Guilloteau et al. (1997). 
They both are clearly detected at 350 /im with flux den- 
sities of 106 ± 7 and 52 ± 8 mJy, respectively. However, 
BRI 0952— 0115, which i s the third z > A quas ar in which 
CO has been measured ( puillotcau et al. 1999 ), is not de- 
tected at 350 fim at a 3 cr level of 65 mJy. This upper limit 
is consistent with the wea k flux density of BRI 0952—0115 
at 1.3 mm 2.8 ± 0.6 mJy ( pmont et al. 1996b| ; puiUotcaij 
et al. 1999), and a temperature of ~ 50 K (see below). The 
radio-quiet quasar HM 0000— 263 at z = 4.11, wh ich was 
not detected at 1.25 mm using the 30m telescope ( pmont 
et al. 199"6b| ), due to its low dechnation, shows a large flux 
density at 350 /zm (134±29mJy). Measurements at other 
wavelengths would be useful to further constrain the prop- 
erties of this object. 

The detection of the z = 3.8 radiogalaxy 4C41.17 with 
a flux density of 37 ± 9 mJy at 350 /Ltm is one of the most 
sensitive measurements of this study. This sensitivity was 
reached after only 3/4 of an hour of on-sourcc time and 
deflnes the limits which can be achieved with SHARC in 
the p ointed observing mode under excellent weather condi 



tions. I A marginal detection (4 cr) was achieved at 350 /^m 
of PC 2047±0123, a. z = 3.80 quasar studied by |Ivisoif " 
(1995) . Finally, the 350 /im flux density of the Clover- 
leaf (H 1413±117) is significantly higher than the value 



publi shed by Barvainis, Antonucci, fc Coleman (1992) 
i.e. 293±14mJy as compared to 189±56mJy. 



226 ± 34 mJy in excellent agreement with the measure- 
ment of 224 ± 38mJy at 438 /itm of Barvainis, Antonucci, 
& Coleman (1992), as shown in Figure ||. 

A greybody was fit to the data points S^, for wave- 
lengths of 350 /.tm and longward, as a function of the rest 
frequency f — Vohs/ (1 + z), of the form 



= BM[l - exp(-T)] with r = {ly/iyoY 



(1) 



where — 2.4 THz is the critical frequency at which the 
source becomes optically thin and H. is the solid angle of 
emission. The shape of the fitt ed greybody is very w eakly 
dependent on the value of (Hughes et al. 1993). The 
data were each weighted by their statistical errors in the 
minimization. This yields the dust temperature, dust mass 
(following Hildebrand (1983), using a dust mass emission 
coefficient at i'q of 1.9 m^ kg""'^), and luminosity of the 
sources. When /3 is considered as a free parameter, we find 
that the average value of /3 = 1.5 ± 0.2 for the detected 
sources. The fits shown in Figure ^ assume an emissivity 
index oi (3 = 1.5. We estimated the 1 a uncertainty in 
the temperature by examining the xt hypersurface in the 
range 1 < /3 < 2, similarly to the method of Hughes et al. 
(1993). To evaluate the uncertainties associated with the 
mass and luminosity, derived from the fitted temperature 
and /? = 1.5, we used the maximum and minimum values 
of the mass and luminosity which are compatible with the 
data plus or minus the statistical error. No lensing am- 
plification was taken into account. The temperature, dust 
mass and luminosities derived under these assumptions are 
given in Table ^ Two of the sources with upper limits have 
1.25mm detections ( pmont et al. 19961:| ), which, together 
with the 350 /im data, yields an upper limit to their tem- 
perature. For Q0842 + 3431, we find that Tdust < 40K 
while for BR 1117 - 1329 a limit of Tdust < 60K is found. 
If the dust is at the temperature limit, these quasars have 
dust masses < lO^M©. For the other sources, an estimate 
of the maximum luminosity has been given under the as- 
sumption that each object has a temperature of 50 K and 
an emissivity index of (3 = 1.5. 

The total luminosity is probably underestimated, since 
a large luminosity contribution from higher temperature 
dust cannot be ruled out for most sources. However, 
in the case of H 1413+117 and IRAS F10214+4724, the 
available IRAS data allow us to fit an additional warm 
component. For IRAS F10214+4724, the cold compo- 
nent model carries roughly 60% of the total luminosity; in 
the case of H 1413+117, which has a hotter mid-IR spec- 
trum, the total luminosity is underestimated by a factor 
of 3. Under the assumption that the majority of the lu- 
minosity is carried by the cold component (Table |l|), the 
median luminosity-to-mass ratio is around 100 Lq/Mq, 
assuming a gas-to-dust ratio of ~ 50 similar to that of 
IRAS F10214+4724 and H 1413+117 (iDowncs et al. 1992 



Barvainis et al. 1995) or of ul traluminous infrared galax- 
ies (ULIRGs), i.e. 540 ± 290 ( [Sanders, ScoviUe, fc Soifa| 



1991). The peak emission in the rest frame is found to 
be in the wavelength range Apeak ~ 60 — 80 /im (Figure |^) 
implying dust temperatures of 40-60 K (Table ll). These 
temperatures are nearly a factor of two lower th an prev i- 



ously estima ted for ultraluminous sources (e.g. phini &| 
We have Kriigel 1994). If the temperature range we find is typical 



also obtained for the Cloverleaf a 450 /Ltm flux density of for the cold component of highly redshifted objects, multi- 
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Fig. 2. — Spectral energy distributions of six of the high-redshift objects discussed in this paper. The 350 /im fluxes, 
shown as squares, are the measurements made with SHARC at the CSO. Fluxes taken from the literature are shown as 
circles - see Tables 1 for references. The observed fluxes at Aobs > 350 ^m were approximated by grey body spectra, shown 
as solid lines, assuming an emissivity index P = 1.5. The dust temperature results from a non-linear least-squares fit to 
the data - see text. 



band photometric studies in the submillimeter/FIR, such 
as planned with FIRST, will provide reasonably accurate 
redshift estimates for the sources detected in deep field 
surveys. 

The global star formation rate in each of the detected 
sources can be estimated using the relation of Thronson 

& Telesco (1986): SFR - lO-i^LpiR/io/iroo ^^0 y^"^ 
with ^ '--^ 0.8-2.1. For our mean luminosity of 1.7 x 
1013/1-2 Lq, this yields a SFR of - 2000 Mq/iJ^q^o yr^^ 
(uncorrected for lensing) if all the submillimeter flux is 
from a starburst component. If we assume a final stel- 
lar mass of 2 X IO^^Mq, a value appropriate to a gi- 



ant elliptical like M87 ( pkazaki fc Inagaki 1984 ), then the 
timescale for formation in a single massive starburst is 
~ lO^ft.^00 yi'- Given the large mass of dust already present 



in these quasars, a substantial amount of this star forma- 
tion must already have occurred. For the most distant 
quasars, where the age of the universe is similar to the 
derived formation timescale, this implies a very high red- 
shift (z > 5) for the era of initial star formation, in agree- 
mcnt with models of hig h redshift Lyman-alpha emitters 
(Haiman & Spaans 1999). 
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